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ABSTRACT 

The number of methods used to study the properties of galaxies is increased, and 
testing these methods is very important. Galactic globular clusters (GCs) provide an 
excellent medium for such test, because they can be considered as simple stellar pop- 
ulations. We present ages and metallicities for 40 Galactic GCs as determined from 
three publicly available techniques, including colour, Lick-index and spectrum-fitting 
methods, based on Bruzual & Chariot evolutionary population synthesis (EPS) mod- 
els. By comparing with the ages obtained from colour-magnitude diagrams (CMDs) 
and metallicities obtained from spectra of stars, we are able to estimate the ability of 
these methods on determination of GCs parameters, which is absolutely necessary. As 
a result, we find that: (i) for the metallicity, our derived metallicities agree with those 
derived from the spectra of stars, Lick-index method is suitable to study metallicity 
for the stellar population systems in the range of —1.5 <[Fc/H]< —0.7 and spectrum- 
fitting method is suitable to study metallicity for the stellar population systems in the 
range of —2.3 <[Fc/H]< —1.5; (ii) for the age, these three methods have difficulties in 
age determination, our derived ages are smaller (about 2.0 Gyr, on average) than the 
results of CMDs for all these three methods. We use Vazdekis and Maraston models to 
analyze whether our results are dependent on EPS models, and find that the tendency 
of these two models is the same as that of Bruzual & Chariot models. Our results are 
independent of the EPS models. In addition, our test is based on the old GCs and our 
conclusions may hold for old stellar population systems. 

Besides the age-metallicity degeneracy, we examine the possible effects of other 
factors (horizontal branch morphology, blue straggler stars, binary interactions and 
a— enhancement) and give quantitative analysis of the influences of these factors on age 
determinations (except for a— enhancement). For colour and spectrum-fitting meth- 
ods, the age can be underestimated about 0.0 — 3.0 Gyr, 0.0 — 2.0 Gyr, and 0.0 — 3.0 Gyr 
due to influences of horizontal branch, blue straggler and binary stars, respectively. 
And for Lick-index method, the lower limit of maximal change of age is 6.0 Gyr, 5.0 Gyr 
and 3.0 Gyr due to influences of horizontal branch, blue straggler and binary stars, 
respectively. 

Key words: Galaxy: globular cluster: general - galaxies: fundamental parameters - 
galaxy: stellar content. 



1 INTRODUCTION 

Understanding the formation and evolution of galaxies is 
one of the important questions of cosmology. With the in- 
creased telescope size and improved instruments, we can 
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obtain multiband photometric and high-resolution spectro- 
scopic data of extragalactic sources. People can take advan- 
tage of these data to study the properties of galaxies to a 
maximum extent. 

Many techniques have been emerged to study the prop- 
erties of galaxies, and evolutionary population synthesis 
(EPS) is one of these techniques which has been devel- 
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oped rapidly since it was first introduced by Tinsley (1968). 
By EPS method, we can compare the observational data 
of clusters and galaxies with EPS models to investigate 
their properties. So far many EPS models (e.g. Vazdekis 
1999; Schulz et al. 2002; Bruzual & Chariot 2003, here- 
after BC03; Le Borgne et al. 2004; Zhang et al. 2005a,b; 
Maraston 2005; Zhang & Li 2006) emerged and were applied 
widely in the study of galaxy properties. Based on these 
EPS models, people adopted different methods to study 
the properties of stellar population (SP) systems, e.g. using 
colours to estimate age and metallicity (Dorman, O'Connell 
& Rood 2003; James et al. 2006); choosing suitable Lick/IDS 
absorbtion-line indices to break the age-metallicity degener- 
acy (Worthey 1994; Worthey & Ottaviani 1997; Trager et al. 
1998) ; using full spectrum-fitting to match SP system spec- 
tra with model spectra (Koleva et al. 2008, hereafter K08) 
and so on. 

Now it becomes the most important to test them and 
estimate the reliability for these methods, especially for the 
fact that a lot of extragalactic studies are heavily dependent 
on these methods. Galaxies are not the ideal test medium for 
these methods because of their complex mixture of SPs. On 
the contrary, Galactic globular clusters (GCs) are suitable 
test beds for this purpose. Presumably formed in a single 
burst, GCs can be characterized by a single age (t) and 
a single metallicity (Z), they are considered as the simple 
stellar populations (SSPs). The age and metallicity can rep- 
resent the overall SP component of a GC. Furthermore, the 
ages and metallicities of them can be known from spatially 
resolved observations (such as colour-magnitude diagrams, 
hereafter CMD and spectroscopy of individual stars), which 
provide us a benchmark for this test. Therefore we use the 
Galactic GCs to study the reliability of these three methods. 

Similar to this work, GCs are also used in other tests, 
such as testing different models and methods based on GCs 
(Mendel, Proctor & Forbes 2007; Riffel et al.2011; K08). 

• Mendel et al. (2007) have used Galactic GCs to test 
three EPS model predictions of age, metallicity and a- 
element abundances. Riffel et al. (2011) have employed the 
near-infrared integrated spectra of 12 Galactic GCs to test 
near-infrared EPS models. 

• Meanwhile, K08 has used Galactic GCs to investigate 
the consistency and reliability of using full spectrum-fitting 
on the determination of SSP equivalent age and metallicity. 

Beyond the work of K08, we use Galactic GCs and BC03 
models to investigate the reliability of three (colour, Lick- 
index and spectrum-fitting) methods on the estimations of 
age and metallicity. Through comparing these results with 
those obtained by other techniques (such as CMDs, spec- 
tra of stars) , it can provide a valuable examination of these 
derivative methods of EPS models, and give a guidance to 
users of these methods. 

The outline of this paper is as follows. We give the 
Galactic GC sample data, population synthesis models and 
stellar spectral library in Section 2. We analyze the methods 
used in this work in Section 3. In Section 4, we present the 
results and discussions of this work. The influences of hori- 
zontal branch (HB) morphology, blue straggler stars (BSs), 
binary interactions and a— enhancement on the age deter- 
minations are given in Section 5. The summary is presented 
in Section 6. 



2 THE GALACTIC GC SAMPLE DATA, 

POPULATION SYNTHESIS MODELS AND 
STELLAR SPECTRAL LIBRARY 

At first, we describe the Galactic GC sample data and the 
EPS models in this Section. In order to discuss the influences 
of HB morphology, BSs and binary stars on the derived ages, 
we introduce the high spatial resolution photometric data 
and stellar spectral library. 



2.1 The sample data 

In this work, we choose 40 Galactic GCs (which have both 
photometric and spectroscopic data) as our sample and use 
the observed data from different studies. The spectroscopic 
data of Galactic GCs are taken from the spectral library of 
Schiavon et al. (2005, hereafter S05). These spectra are ob- 
tained from the Cerro Tololo Inter- American Observatory 
Blanco 4-m telescope with the Ritchey-Chretien spectro- 
graph. And these spectra cover the wavelength 3350-6430 A 
with a resolution of about 3.1 A. 

The photometric data of our Galactic 
GC sample are taken from the 2010 ver- 
sion of the Harris' catalogue (Harris 1996; 
http:/ /physwww.physics.mcmaster.ca/~harris/mwgc.dat ). 
The selected targets span a wide range of cluster pa- 
rameters, including metallicity, HB morphology, Galactic 
coordinates and Galactocentric distance. And the main 
characters of them are listed in Table Al, which are taken 
from Harris' catalogue. 

It is well known that the ages of GCs can be obtained 
with much higher accuracy by CMDs |De Angeli et al.l 
2005), and the metallicities can be obtained by the high 
resolution spectra of stars l|Carretta et al.l l2009h . In or- 
der to give a comparison, w e take the ages of 30 GCs 
from Salaris fc Wei s s (20 02), |Pe Angeli etal] (|2005l ) and 
iMarin-Feanch et all (|2009l '). which determine ages of GCs 
from the CMDs, an d the ages are mainly taken from 
ISalaris fc Weiss! ( 
are taken from 



20021). The m e tallici ti es of all these 40 GCs 
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(|l997f) . iKraft fc Ivans! l|2003l ) and ICarretta et all (|20QgF T 



which determine metallicities of GCs from sp ectra of stars, 
and the metallicities are mainly taken from IZinn fc West! 
l|l984h . 



2.2 The stellar population synthesis models 

To test these three methods on the parameter determi- 
nations, we adopt the standard BC03 models. The stan- 
dard BC03 population synthesis models utilize the Padova 
1994 evolutionary tracks, the STELIB library (Le Borgne 
et al. 2003) and the Chabrier (2003) initial mass func- 
tion (IMF) with stellar mass limits of 0.1 and 100 M Q . 
This set of models presents the intermediate-resolution 
spectra (~3A), some spectral feature indices and colours 
(based on various systems, including Johnson- Cousins and 
AB systems). The SSP models given by BC03 span a 
wide range of age r (0.1 Myr — 20.0 Gyr) and metallicity 
Z (0.0001—0.05), but these models provide six metallici- 
ties: [Fe/H] = -2.3, -1.7, -0.7, -0.4, +0.0, and +0.4 (here 
we adopt Z Q =0.02). 

In order to discuss the influence of binary stars on age 
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determinations in Section 5.3, we also adopt the other two 
EPS models, one is the models of SSPs (Zhang et al. 2004) 
and the other is the models of binary stellar populations 
(BSPs, Zhang et al. 2005a). These models were b uilt on the 
basis of the Cam bridge stellar evolution tracks l|Eggleton I 
Il97ll . 1 19721 . I1973T ). BaSeL-2.0 stellar atmosphere models 
(Lejeune, Cuisinier & Buser 1997, 1998) and various ini- 
tial distributions of stars. The main input parameters of 
the BSP models are as follows: (i) the initial mass of the 
primary is chosen fr om the app roximation to the IM F of 
iMiller fc Scalol (|l979h as given bv lEggleton et alj (| 19891 ): (ii) 
the initial secondary-mass distribution, which is assumed to 
be correlated with the initial primary-mass distribution, sat- 
isfies a uniform value, i.e., n(q) = l; (iii) the distribution of 
orbital separations is taken as constant log a for wide bi- 
naries and falls off smoothly at close separations; (iv) the 
eccentricity distribution satisfies a uniform form: e = X. 
Other assumptions and more details see Zhang et al. (2004) 
for SSP models and Zhang et al. (2005a) for BSP models. For 
those previous m odels, Zhang e t al. adopt t he 2.2-version of 
BaSeL library of lLeieune etaD l| 19971 . 1 19981 ). which is a low- 
resolution library. However, the spectra of S05 library are 
intermediate resolution and we calculate the Lick indices by 
degrading to the resolution of Lick system (see Section 4.1), 
we do some modification for these two models, instead of us- 
ing BaSeL library we choose the high-resolution (~0.lA, 
see Section 2.3) Bluered library of Bertone et al. (2008) in 
this work. And other input physics remain unchanged. 

2.3 Other photomethic data and stellar spectral 
library 

For the purpose of investigating the influences of HB mor- 
phology and BSs on age determination, in the first step we 
should select the HB stars and BSs f rom CMDs of GCs. 
CMDs taken from iPiotto et all l|2002j . hereafter P02) are 
used to obtain HB stars and BSs for GCs. These data are 
high spatial resolution photometric data and observed with 
the Wide-Field Planetary Camera 2 on board of Hubble 
Space Telescope (HST/WFPC2) camera in the F439W and 
F555W bands. The HST snapshot catalogue contains pho- 
tometric tables for 74 GCs and the B and V magnitudes are 
also obtained from the HST photometry through the trans- 
formations given by P02. There are 26 clusters in common 
between the spectral library of S05 and P02 photometric 
catalogue. 

In the second step we should give the corresponding 
colours and spectra for each HB star and BS. In order to 
do this, we adopt two stellar spectral libraries. One is the 
2.2-version of BaSeL library of Lejeune et al. (1997, 1998), 
which is a low-resolution library. This library provides ex- 
tensive and homogeneous grids of low-resolution theoret- 
ical flux distributions in the range of 91-1600 000 A and 
synthetic UBVRIJHKLM colours for a large range of stel- 
lar parameters: 50 000 > T cff /K ^ 2000, 5.50 > logg ^ -1.02 
and 1.0^ [Fe/H] ^ —5.0. By this library we can obtain the 
colours for HB stars and BSs. In order to obtain high- 
resolution spectra for them we adopt the other library, 
it is Bluered library of Bertone et al. (2008), which is 
a high-resolution (R > 500 000) library. This library con- 
tains 832 theoretical stellar spectra covering the optical 
range of 3500— 7500 A (see Bertone et al. 2008 for more 



details). It spans a large range of T c s (4000-50,000 K) 
and logg (0.0—5.0) at six metallicity values ([Fe/H] = 
-3.0, -2.0, -1.0, -0.3, +0.0 and +0.3). 



3 METHODS 

Having described the sample data and EPS models, we now 
begin to introduce the three methods used in this paper, 
including colour, Lick-index and spectrum-fitting methods. 

3.1 Colour method 

As shown by Yi et al. (2004), the U - B and B — V colours 
can be well used to derive the ages and metallicities for 
GCs. Anders et al. (2004) also analyzed the reliability and 
limitation of the combination of various passbands and they 
found that the U and B bands are important for the age 
and metallicity determinations, and the V band and near- 
infrared data can provide additional constraints. They listed 
the preferable passband combinations in their Tables 1 and 
2, and the UBVI are the preferable combinations for the SP 
systems with age larger than a few Gyr from their Table 
2. So we select three colours (U — B, B — V and V — /) 
from Harris' catalogue to fit with models to obtain the ages 
and metallicities of GCs. But for the other SP systems, the 
selecti on of colours may be different, such as lMaraston et ail 
(2001) have found that the ages of young star clusters can 
be well determined based on B — V and V — K colours, 
and wider photometric data can be use d well to constrai n 
the parameters for distant SP systems l|Barro et al.ll201lT ). 
In this procedure, we set the ages and metallicities of SPs 
free and select the best-fitting model by minimizing the x 2 
following the formula: 

2 1 \ ~> ( Cm.,.™ — Co,„\ 

where the Cm,i,n is the n-th colour corresponding to the i-th 
SSP model. The Co,n and ao,n are the n-th colour and error 
of the Galactic GCs. N^oi is the degree of freedom. Error 
cro,n is not given in the Harris' catalogue, so we assume that 
the error of U — B, B — V and V — I is 0.04 mag. 

Errors on the age and metallicity estimation are deter- 
mined from the x 2 contours. Because we use three colours in 
this work, the A^dof is 2. And for A^of = 2, the la, 2a and 3a 
errors correspond to the formal confidence intervals of 68.7% 
(Ax 2 = 2.3), 95.4% (Ax 2 = 6.17) and 99.73% (Ax 2 = 11.8; 
Avni 1976), respectively. 

Fig. 1 shows an example of the method of using contour 
map of x 2 to obtain age and metallicity for NGC 6235. The 
star stands for the Xmin an d corresponds to the obtaining 
r and [Fe/H]. The shaded area maps la confidence region, 
and the At and A [Fe/H] are obtained based on this shaded 
area. 

3.2 Lick-index method 

Lick/IDS indices have been used to break the age-metallicity 
degeneracy (Trager et al. 2000; Terlevich & Forbes 2002; 
Gallazzi et al. 2005) of SP systems. In this paper, 
we select the Balmer (H/3, H8a, H5f, H7A and H7F), 
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Figure 1. Contour map of \ 2 as a function of age r and metal- 
licity [Fe/H] for an example of NGC 6235. From the center to the 
edge the value of \ 2 increases from 1.1 to 3.3 in interval of 0.22. 
The black star stands for the minimum \ 2 an d corresponds to 
the obtaining r and [Fc/H]. We list the parameters of NGC 6235 
on the upper right corner. The gray shadow represent la error 
contour, and the errors are obtained based on this gray shadow. 
The dashed lines mark the BC03 model grids, we do interpolation 
for the mctallicity (see Section 4.1). 




02460 2 40 2 4 

Mgb <Fe> [MgFe] 



Figure 2. H/3 versus three indices (Mg6, <Fe> and [MgFe] ) for 
the models of Thomas, Maraston & Johansson (2011a), the ages 
are from 1 to 15 Gyr and [Fe/H] are from —2.3 to 0.67 that are 
labeled in the figure. The gray solid and red dashed lines stand 
for [a/Fe]= 0.0 and 0.5, respectively. 



Mg 6, Fe5270 and Fe5335 indices to study the parame- 
ters of GCs (Balmer indices are age-sensitive indices, Mg b, 
Fe5270 and Fe5335 are metal-sensitive indices). The H<5 
and H7 are high order Balmer indices which measured 
with a narrower (~ 20 A; marked by subscript F) and a 
wider (~ 40 A; marked by subscript A) central bandpass 
ijWorthev fc Ottaviani|[l997l ). Among these five Balmer in- 
dices, H/3 and H7F are little sensitive to a/Fe and the 
other three indices are sensitive to a/Fe. And in Sec- 
tion 5.4, we will discuss the influence of a-enhancement 
on age determination for Lick-index method. For the 
three metal-sensitive indices, we adopt the following def- 
inition: [MgFe]' = [Mg 6(0.72 x Fe5270 + 0.28 x Fe5335)] 1/2 , 
which defined by Thomas, Maraston & Bender (2003). The 
[MgFe] has been found to be good tracer of metallicity. 
Adopting the models of Thomas, Maraston & Johansson 



(2011a) we present this character in Fig. 2, the grids of 
H/3 versus Mgfe, <Fe> (<Fe> = l/2(Fe5237+Fe5335)) and 
[MgFe]' with two types of a/Fe ([a/Fe]= 0.0 and 0.5), the 
gray solid and red dashed lines stand for [a/Fe]= 0.0 and 
0.5, respectively. From them we can see that the influence 
of a/Fe ratio is obvious for H/3 versus Mgfe and <Fe> grids, 
and is small for H/3 versus [MgFe] grid. The Mgb increases 
and <Fe> decreases with increasing a/Fe, and the combi- 
nation index [MgFe]' counteract the effect of a/Fe ratio on 
Mgfe, Fe5237 and Fe5335. This certify that the [MgFe]' is 
independent o f a/Fe and can be c onsidered as good tracer 
of metallicity dThomas et alj|2003f ). 

In Fig. 3, we show the distribution of GCs on five 
Balmer (H/3, FISa, H5f, H7A and H7F) indices versus 
[MgFe] planes. We can see that the majority of GCs lie 
beyond the grid covered by H/3 (this problem does not in- 
volve all Balmer indices, Thomas et al. 2011a), which would 
be inconsistent with the typical universe age 13.6 Gyr. This 
is a general problem (see Fig. 2 of Mendel et al. 2007) for 
all studies of using H/3 index to study the GCs based on ex- 
isting models. This problem is also clearly demonstrated by 
the 'zero-point' offset (some GCs lie beyond the grid of H/3 
and most GC ages are greater than 14 Gyr). Although these 
problems have been pointed by some groups (Vazdekis et 
al. 2001; Schiavon et al. 2002; Cenarro et al. 2008), but few 
works have been done on studying the specific sources of this 
problem or quantified their effects on the derived ages and 
metallicities. Similar to the technique described in Puzia et 
al. (2005), the ages and metallicities for individual GCs are 
computed as the weighted mean of the parameters derived 
from five Balmer indices versus [MgFe] grids. 

Due to the influence of hot blue HB stars, the isochrones 
for old (r > 8.0 Gyr) stellar populations with metallicity be- 
low [Fe/H] = —0.7 tend to overlap (Maraston & Thomas 
2000; Mendel et al. 2007), this nature is reflected by H/3 in 
panel (a) of Fig. 3. This brings an ambiguity in determining 
ages and metallicities of GCs and unnaturally extends the 
age distributions for old ages for Lick-index method. Note 
that, for the GCs lying beyond the Lick indices grids, we 
assume the values of GCs are 15.0 Gyr. 

3.3 spectrum-fitting method 

With the increase of spectroscopic data, people use spec- 
tra to study the properties of SP systems. So far there are 
many kinds of full spectrum-fitting techniques used in such 
study (e.g. Cid Fernandes et al. 2005; Mathis, Chariot & 
Brinchmann 2006; Tojeiro et al. 2007; K08). 

We employ a modified version of the starlight code by 
Cid Fernandes et al. (2005) to analyze the observed spec- 
tra of Galactic GCs. The starlight code is originally used 
to study the properties of galaxy, and is achieved by fit- 
ting the observed spectrum Fo with a model spectrum Fm 
that mixed by /V* SSPs with different ages and metallicities 
from the BC03 models. The code is carried out with a sim- 
ulated annealing plus Metropolis scheme (Cid Fernandes et 
al. 2001), which searches for the minimum 

X =^[(Fo-fMK] 2 , ( 2 ) 

A 

where cj^ 1 is the error in Fq. The line-of-sight stellar mo- 
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Figure 3. The distribution of GCs on the Balmer indices and [MgFe] planes. Model grides from BC03 are shown: age is constant 
for solid line (top to bottom, 1, 2, 4, 5, 8, 10, 12 and 15 Gyr), and metallicity is constant for dashed line (left to right: [Fe/H] = 
—2.3, —1.7, —0.7, —0.4, 0.0 and +0.4). Line strengths of Galactic GCs are plotted on SSP grides. The red open squares, blue open circles 
and purple open pentacles stand for GCs with red, blue and unknown HBR, respectively. In the top left plane (H/3 versus [MgFe] plane), 
we give these values of age and metallicity. 



tions are modeled by a Gaussian distribution centered at 
velocity i;* and with dispersion <r*. The output includes the 
distributions of stellar age, metallicity, extinction, velocity 
dispersion and stellar mass. 



4 RESULTS AND DISCUSSIONS 

We use colour, Lick-index and spectrum-fitting methods 
based on BC03 model to obtain the ages and metallicities 
for Galactic GCs, and compare them with those determined 
from spatially resolved observations (such as CMDs, spectra 
of stars) . We also compare the ages and metallicities derived 
from other two models (Vazdekis and Maraston models) to 
test whether our results are dependent on EPS models. Be- 
fore displaying the results of this work, we should introduce 
some pretreatments of model, sample data and methods. 
The details are given as follows. 



4.1 Pretreatments 

• For the EPS model, as said above, the BC03 model pro- 
vides six metallicities, we linearly interpolate colours, Lick 
indices and spectra of SSPs that span [Fe/H] from —2.3 to 
+0.4 in increments of 0.1 dex. 

• In order to discuss the influence of HB morphology 
on age determinations in Section 5.1, we divide the GCs 
into three groups according to the value of HB ratio (HBR 



Q) from the Harris' catalogue. One is the blue HBR type 
(HBR >0), the second is red HBR type ( HBR <0) and the 
third is unknown HBR (HBR is not given). 

• For the colour method, we do extinction corrections 
for obser ved magnitude s of G Cs by adopting the extinction 
curves o f ISchlegel etafl (|l99Sh and the E(B - V) of Harris' 
catalogue (the 8th column of Table Al). 

• For the Lick-index method, in this work we calculate the 
Lick indices from the spectra of Galactic GCs by degrading 
to wave-dependent resolution of the Lick system (see Section 
4.4 of BC03), without adopting the results by using fitting 
functions. 

• For the spectrum- fitting method, we adopt a modified 
version of the starlight code. Because it is generally ac- 
cepted that the GC can be represented by a SSP, only one 
component in the base. We fit the observed GC spectrum 
Fo with each of the SSP spectrum Fi, do not fit the ob- 
served spectrum wi th a linear combination of N* SSPs. In 
this study, following lCid F crnandes~ et al.l (|2005l ). we find the 
best fitting SSP, which matches a given observed spectrum 
of GC, through a standard \ 2 minimization procedure: 

A 

instead of equation (2), where Fi corresponds to the spec- 
trum of i-th SSP model. And in this work we construct a 

1 HBR= (B - R)/(B + R + V), B is the number of blue HB 
stars, R is the number of red HB stars and V is the number of 
RR Lyrae variable. 
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Figure 4. Top panels: the comparisons between the metallicities obtained by three methods and the spectra of stars for Galactic GCs. 
Symbols have the same meanings in Fig. 3. Bottom panels are the comparisons between the ages obtained by three methods and the 
CMDs. Left, middle and right panels are comparison results for methods of colour, Lick-index and spectrum-fitting, respectively. 



base of 18 ages x 24 metallicities = 432 SSP model spectra. 
The use of starlight to study the integrated spectra of clus- 
ters has also been extensively discussed by Cid Fernandes & 
Gonzalez Delgado (2010). 

In the spectrum-fitting process, we use the spectra in the 
3700—5700 A range because the red and blue parts of the 
observations have poor quality. Errors on the age and metal- 
licity estimation are also determined by the \ 2 contours. 



4.2 Comparison of parameters with other studies 

Using the methods and pretreatments outlined in Sections 3 
and 4.1, we obtain ages, metallicities and corresponding er- 
rors for GCs. In order to distinguish the ages and metallici- 
ties derived from different methods, we define some param- 
eters in this paper, rc and [Fe/H]c stand for the age and 
metallicity determined from colour method, ri and [Fe/H]i 
represent those obtained from Lick-index method, and rs 
and [Fe/H]s stand for those determined from spectrum- 
fitting method. The parameters of literatures are repre- 
sented by tl and [Fe/H] L . In Table A2 we list the GC IDs in 
the first column. In the 2nd to 7th columns, rc, [Fe/H]c, Tt, 
[Fe/H]i, rs and [Fe/H]s are given with errors, respectively. 
And in 8th and 9th columns we list the rL and [Fe/H]L. 

In the top panels of Fig. 4, we give the comparison be- 
tween our GC metallicities and those given in the literatures 
that obtain metallicities from the spectra of stars. And (a), 
(b) and (c) panels are the comparison results for colour, 
Lick-index and spectrum-fitting methods, respectively. The 
vertical bars are the errors in metallicities derived from three 
methods. From them we see that the metallicities obtained 
by three methods are in agreement with literature value for 



the whole sample. Meanwhile we also find that all GCs with 
blue HBR have lower metallicity than those with red HBR 
for these three methods, which agrees with the study of Lee, 
Yoon & Lee (2000). Our specific analysis are given as fol- 
lows. 

(a) For colour method, [Fe/H]c agrees with [Fe/H]L for 
the whole sample, but there exists discreteness. We can see 
that [Fe/H]c is a bit smaller than [Fe/H]L in the range of 
[Fe/H] L < -1.0 and [Fe/H]c agrees with [Fe/H] L for the 
range of [Fe/H] L > -1.0. But for the value of [Fe/H] L is 
about -1.1, there exist three GCs (NGC 6171, 6342 and 
6652) having lager [Fe/H] c than [Fe/H] L . From the [Fe/H] c 
and rc of these GCs listed in Table A2, we see that rc 
is smaller than rL, so these GCs can be affected by age- 
met allicity degeneracy. 

(b) For the Lick-index method, the [Fe/H]i is perfectly in 
agreement with [Fe/H]L in the range of —1.5 < [Fe/H] < 
—0.7. But there have difficulties in metallicity determination 
for the range of [Fe/H] < —1.5, because it tends to overlap 
for Balmer and [MgFe] indices at low-metallicity range for 
BC03 model. And we can see this phenomenon clearly in the 
panel (a) of Fig. 3. This indicates that Lick-index method is 
suitable to study metallicity for SP systems in the range of 

— 1.5 < [Fe/H] < —0.7 and has some difficulties in studying 
SP systems for the range of [Fe/H] < —1.5. 

(c) For the spectrum-fitting method, we find that [Fe/H]s 
matches [Fe/H] L in the range of -2.3 < [Fe/H] < -1.5, and 
[Fe/H] s is smaller than [Fe/H] L for the range of [Fe/H] > 

— 1.5. All these results imply that the spectrum-fitting 
method may be suitable to study metallicity for SP sys- 
tems in the range of —2.3 < [Fe/H] < —1.5 and [Fe/H]s may 
be smaller than [Fe/H]L for the range of [Fe/H] > —1.5. 
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On the whole, our metallicities obtained from three 
methods match those determined from spectra of stars in the 
entire metallicity range spanned by the GCs. The Lick-index 
method is suitable to study the metallicity in the range of 
— 1.5 < [Fe/H] < —0.7 and spectrum-fitting method is suit- 
able to study the metallicity in the range of —2.3 < [Fe/H] < 
-1.5. 

In the bottom panels of Fig. 4, we compare our derived 
GC ages of three methods with the values of literature that 
determine ages from CMDs. And (d), (e) and (f) panels 
are corresponding to colour, Lick-index and spectrum-fitting 
methods, respectively. From them we can see there exists 
large discrepancy in age determinations, and the errors are 
also large for these three methods. Our specific analysis are 
given as follows. 

(d) For the colour method, most of the GCs with blue and 
unknown HBR have relative lower rc than tl, but there has 
no obvious tendency for those with red HBR. The age of 
GCs with blue HBR can be influenced by HB stars for the 
existence of HB stars making GCs look younger. Five GCs 
are found to possess extreme low rc, including two with un- 
known HBR (NGC 6388 and 6441) and three with red HBR 
(NGC 6171, 6342, and 6652). Among these GCs, three have 
relative lager [Fe/H]c than [Fe/H]L, which can be affected 
by the age- metallicity degeneracy (see Fig. 4a). Meanwhile 
from the 4th and 5th columns of Table Al, we know that 
the other two GCs are close to the Galactic bulge, which are 
badly contaminated by the field stars or affected by differ- 
ential reddening, so we can not get reliable parameters for 
them. 

(e) For the Lick-index method, we can see there exist some 
GCs with extreme large ri (about 15.0 Gyr), and except 
these GCs most of GCs have smaller ri which are influenced 
by HB stars, especially for those with blue HBR. About 
eight GCs have large ri of 15.0 Gyr which lie beyond the 
grids, including two with red HBR (NGC 1851 and 6637), 
one with unknown HBR (NGC 5946) and five with blue 
HBR (NGC 1904, 2298, 3201, 6254 and 7078). Similar to 
above, Lick-index method has some difficulties in parameter 
determinations for the range of [Fe/H] < —1.5, and most of 
these eight GCs are metal-poor ([Fe/H]i< —1.5), so we ob- 
tain extreme large n for them. 

(f) For the spectrum-fitting method, most GCs with blue 
HBR have relative smaller rs than tl, but there has no ob- 
vious tendency for those with red and unknown HBR. From 
panel (c), we know that [Fe/H]s is consistent with [Fe/H]L 
in the range of [Fe/H] < —1.5 and most of GCs with blue 
HBR are metal-poor, so the rs of GCs with blue HBR are 
affected by HB stars rather than age- metallicity degener- 
acy. For the GCs with red and unknown HBR, one part has 
large rs which can be influenced by age-metallicity degener- 
acy, and the other part has small rs which may be influenced 
by HB stars and BSs. Note that, relevant to this test is that 
CMD-derived ages also can carry their own problems for this 
test, because the CMD-derived ages depend on the adopted 
tracks and on whether element rati os be taken into account, 
which has also been discussed in iMaraston fc Stromback 
|201ll ). 

On the whole, all of these three methods have difficul- 
ties in age determination and we can not directly say which 
method is better on age determination. Except the influ- 
ences of CMD-derived ages and age-metallicity degeneracy 



in the entire age range, there exist some uncertainties in 
age determinations, some factors (e.g. HB stars, BSs, binary 
stars and a-enhancement) can affect the age determinations 
and we will discuss them in the next Section. 



4.3 Consistency checks 

We u se Vazdekifl (the update d v 9.1 ve rsion, IVazdekis et al.l 



l20ld: 



Falcon-Barroso et al 



EoTTI) 



and 



Maraston 
investigate 



( Maraston fc Stromback! 1201 lh models to 
whether our results are dependent on EPS models. 
The Vazdekis models use t he P adova 2000 theoreti- 
cal isochrones dGirardi et alj bood ), the MILES library 
ijSanchez-Blazquez et al. 200d l and four types of IMF 
with stellar mass limits (Vazdekis et al. 2003) of 0.1 
and 100 Mq. This version covers the wavelength range 
of 3540 -7410 A with a spectral re solution (FWHM) of 
~2.5A (|Falc6n-Barroso et al.l l201ll ). seven metallicities 
0.0001 ^ Z ^ 0.03 and 50 ages across the range of 
0.0630— 17.180 Gyr. In this w ork, we select the unimodal 
Salpeter IMF (|Salpeterl \wEEi) to study the parameters of 
GCs for spectrum-fitting and Lick-index methods. For the 
colour method, because the MILES library has relatively 
limited wavelength coverage and does not extend into 
I-band, we select the c olours from previous models of 
Vazdekis (|Vazdekislll999T ) to study the parameters of GCs. 
This does not influence our results (private discussion). 

The Maraston models assume a stellar evolution 
prescription, which consists of the isochrones and stel- 
lar tracks by Cassisi, Castellani & Castell ani (1997) for 
ages larger than ~30.0Myr and by Ganeva |Schaller et al.l 
1992) for younger populations. Sets of models have been 
also computed with P adova stellar evolutionary mod- 
els (|Girardi et alj kOQO). Four different libraries of flux- 
calibrated empirical stellar spectra and three types of 
IMF have bee n conside r ed. W e select the MILES li- 
brary and the IChabrierl (|2003l l IMF to study the pa- 
rameters of GCs. These updated models cover the wave- 
length range of 350 0— 7430 A with a s pectral resolution 
(FWHM) of ~2.54 A (|Beifiori et al.ll201ll ) and five metallic- 
ities 0.0001 ^ Z 0.04 with different ages across the range 
of 0.0060-15.0000 Gyr (Maraston & Stromback 2011, and 
references therein). Just as said above, w e use the updated 
models of IMaraston fc Stromback! l^Oll! ) based on MILES 
library for spectrum-fitting and Lick-inde x methods. And fo r 
the colour method we select the models of M araston I (|2005h . 
which based on the BaSeL library, to study the parameters 
of GCs. In Table 1 we list the comparison of main model 
ingredients in each EPS model. The first, second and third 
rows list the stellar evolution track, stellar spectral library 
and IMF of the EPS models that we used in this study. The 
4th and 5th rows show the metallicity (Z) and age range 
that these models cover, and the number of metallicities is 
also given in the parenthesis of the 4th row. 

We use these two EPS models to obtain the parame- 
ters of GCs based on three methods described above and 
compare the results of these two models with that of the 
literatures. In Fig. 5, we show the comparisons between the 
parameters of GCs obtained from tree methods based on 

2 http://miles.iac.es/ 
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Table 1. Comparison of main model ingredients in each EPS model. The first, second and third rows list the stellar evolution track, 
stellar spectral library and IMF of the EPS models. The 4th and 5th rows show the metallicity (Z) and age range that these models 
cover, and the number of metallicities is also given in the parenthesis of the 4th row. 



Models 



BC03 



Vazdekis 



Maraston 



Stellar evolution track 
Stellar spectral library 
IMF 

Z (number) 
Age(Gyr) 



Padova 1994 

STELIB 

Chabrier 

0.0001 - 0.0500(6) 
0.0001 - 20.0000 



Padova 2000 

MILES 

Salpeter 

0.0001 - 0.0300(7) 
0.0630 - 17.1800 



Cassisi+Geneva 
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Chabrier 

0.0001 - 0.0400(5) 
0.0060 - 15.0000 
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Figure 5. Comparisons between the results from three EPS models and the literatures. Top and bottom panels are the comparisons 
for metallicity and age, respectively. The open squares, pentacles and circles stand for the results of BC03, Vazdekis and Maraston 
EPS models, respectively. Left, middle and right panels are comparison results for methods of colour, Lick-index and spectrum-fitting, 
respectively. The dispersion (rms) of each model is also listed in the lower right corner of each panel. 



three EPS models (BC03, Vazdekis and Maraston) and the 
literatures, in which the open squares, pentacles and circles 
stand for the results of BC03, Vazdekis and Maraston EPS 
models, respectively. The dispersion (rms) is also listed in 
the lower right corner of each panel, rmsb C 03, rms V azd and 
rmsmara stand for the rms of BC03, Vazdekis and Maras- 
ton EPS models, respectively. We find that the metallici- 
ties obtained from three methods of these two EPS models 
(Vazdekis and Maraston) have an agreement with those of 
the literatures in the entire metallicity range. From the top 
panels (a, b and c panels), we find that the comparisons be- 
tween metallicities obtained by three methods and the spec- 
tra of stars for these two models are similar to that of BC03 
models. For Lick-index method in panel (b), we see that the 
rms vaz d and rms mara are smaller than rms bc 03, this may be 



due to the adoption of MI LES for these two models an d 
STELIB for BC03 models jMaraston fc Strombackll201ll ). 
For the spectrum- fitting method in panel (c), the [Fe/H]s of 
Vazdekis model is larger than that of the other two models, 
so the rmSvazd is smaller than rmsbece and rms mara . From 
these rms, we can know that the dispersion of colour method 
is larger than that of other two methods. All these can not 
influence our conclusions, Lick-index method is suitable to 
study metallicity in the range of —1.5 < [Fe/H] < —0.7 and 
spectrum-fitting method is suitable to study metallicity in 
the range of —2.3 < [Fe/H] < —1.5. From bottom panels, 
we can see that all these three methods have difficulties 
in age determinations for these two EPS models and the 
rms is larger than that of metallicity. There also exist some 
GCs with extreme large t\ (about 15.0 Gyr) for these two 
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Figure 6. Taking NGC 7078 as an example to illustrate the 
method of selecting HB stars and BSs from the CMD. The HB 
stars are marked with red asterisks and BSs are marked with blue 
circles. 

models, and most of GCs have relative smaller re, tj and 
rs than tl, which is similar to that of BC03 models. The 
comparisons of derived metallicities and ages of individual 
GCs show that these three models make different parame- 
ter predictions. However, the whole tendency for these three 
methods of these three EPS models is nearly the same. This 
indicates that our conclusions are independent of the EPS 
models. 



5 UNCERTAINTIES IN AGE 

From the above analyses, we know that our metallicities 
are consistent with the values of literatures within the er- 
ror bars, but there exist some discrepancies in age deter- 
minations. It is well known that the HB stars, BSs, binary 
stars and a— enhancement exist in GCs, however they are 
not considered in most of the EPS models and they are 
important for models. All these can affect the parameter 
determinations for GCs. So in this Section, we analyze the 
effects of some factors (HB morphology, BSs, binary star 
and a— enhancement) on age determinations. 

5.1 Horizontal Branch morphology 

It is widely known that HB morphology can affect the age 
determination (de Freitas Pacheco & Barbuy 1995; Lee et al. 
2000; Maraston & Thomas 2000). The HB stars (especially 
the blue HB stars) can be hot and their existence can mimic 
young or intermediate age populations in old SP systems. 
Hot extended HB stars have been observed in GCs with 
HST by Piotto et al. (1999). 

The main goal of this part is to analyze the contribu- 
tions of observed HB stars to the age determinations for 
three methods. At first, we should select HB stars from the 
CMD and then give corresponding colours, spectra and Lick 
indices for HB stars. We use the distribution of stars in 
the CMD of P02 and select 18 GCs (include five with blue 
HBR, four with unknown HBR and nine with red HBR), 
which have in common with our sample and show obvious 



HB stars, to study the influence of HB stars on age de- 
terminations. We select HB stars according to the method 
described in Gratton et al. (2010), and take NGC 7078 as 
an example to show the selection of HB stars in Fig. 6 (red 
asterisks are HB stars). 

The total observed quantity (Xtot) of a GC includes 
two components: the total HB stars component (Xhb) and 
the component pf(tot-HB)) of other stars in the GC. After 
selecting of HB stars and analyzing Xtot and ^(tot-HB)i we 
can study the influence of HB stars on age determinations. 
The observed spectra and colours of these two components 
(Xhb and X( to t-HB)) can be constructed by follows. 

(i) The integrated colours of these two components, taking 
(B — V)hb and (B — V)( to t-HB) as examples, can be ex- 
pressed as 

(B~V) W = -2.5 log 10 -o.^, HB ^d (4) 

10 -0.4B_y^ w HB 10 -0.4i> iiHB 

(B - VQ ( t t- HB ) = -2-5 log io _ . 4y _gg B 10 - . 4%HB , (5) 

where Nub is the number of HB stars, 6?,hb and v ;,hb are 
the magnitudes for the i-th HB star which given by P02, and 
B and V are the magnitudes for each GC which obtained 
from Harris' catalogue. 

(ii) The spectra of F H b and Fftot-HB) can be obtained by 
following 

«HB - 

ft B = £ / HB and (6) 

t=i 

-F(tot-HB) = -ftot — -Phb , (7) 

where Ftat is the spectrum taken from S05, and / HB is the 
spectrum of i-th HB star which is obtained from the Blubred 
library. In this procedure, we transform the surface flux (/ s l ur ) 
from Bluered library into observed flux (/ HB ) for each HB 
star by adopting the distance of GCs from the sun by Harris 
catalogue. 

For the colours and spectra of HB stars, we can derive 
them from the stellar spectral library based on HB star s 
[Fe/H], T c ff and logg. Because the GC is taken as a SSP, 
the same metallicity (i.e. [Fe/H]L in Table A2) of GC is 
adopted for HB stars for a certain GC. The parameters of 
T e ff and \ogg are obtained as follows. 

(i) The T e g can be obtained from the theoretical 
BaSeL library by the bv magnitudes of P02 and metallici- 
ties of GCs. From the theretical colour— T e g calibration of 
BaSeL library, we can assign the T e g for each HB star. 

(ii) The logg is not easy to obtain, in this work we select 
a bimodal logg for HB stars, a median value of logg = 2.5 for 
HB stars with T e ff<8000 K (include all red and part blue 
HB stars) and logg = 4.0 for HB stars with T efI >8000 K 
(consist of blue HB stars). The reasons are given as follows. 

(a) Moni Bidin et al. (2007) studied the T cff , logs, he- 
lium abundances and masses for HB stars in NGC 6752. And 
they found that the logg of all HB stars covered t he range of 
2.5 < logg < 5.7. Fig. 3 of iMaraston et al. I (|2003l ) displayed 
the theoretical isoc hrones, and the logs range for HB st ars 
was about 2.0 - 4.0. iDorman I fll992h and lLei et all (|201ll . in 
preparation) have presented a large grid of HB evolution se- 
quences (including the evolution of extreme HB stars), and 
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Figure 7. Theoretical Balmer index— T c g relations obtained from the Bluerbd library for different log g (solid, dashed, dotted, 
and dash-dot-dot-dot lines stand for log g— 3.0,4.0,4.5 and 5.0, respectively) at [Fe/H]= —1.0. 



from their data we find the logg range for red HB star is 
about 1.8 — 3.0 and for blue HB star is about 3.0 — 5.5. 

(b) The Balmer indices have a strong dependence on 
the logg for T e g >8000 K, they show decrease with increas- 
ing T e ft within a fixed metallicity. This is shown in Fig. 7, 
from it we see that the differences are large for different logg 
with T cff >8000 K. So for the HB stars with T cff <8000 K, 
the adoption of logg has little effect on Balmer indices, but 
for those with T e ff >8000 K, the adopti on of logq should be 
reason able. In addition, from Fig. 1 of iRecio-Blabco et al.l 
(2006), we can see that all red and part blue HB stars 
are with T e ff<8000 K and part blue HB stars are with 
r e ff^8000 K. Combined with logg range of red and blue HB 
stars, we choose a bimodal logg for them and the boundary is 
T e ff ~8000 K, \ogg = 2.5 (the median value of red HB stars) 
for those with T^ff <8000 K and \ogg = 4.0 (the median value 
of blue HB stars) for those with T cff >8000 K. 

Based on [Fe/H], T c a and logg determined as described 
above and the stellar spectral library, we can obtain the Fhb 
and F( tot _ H B)- Using the methods described in Section 3 we 
can study the influence of HB stars on age determinations. 

For the colour and spectrum-fitting methods, we show 
the influences of HB stars on age determinations in Fig. 8. 
For the sake of simplicity, r represents the age of GC with 
HB stars (correspond to re and rs in Section 4.2) and 

TVtot hb) represents the age of GC without HB stars. We 

fix [Fe/H]c and [Fe/H]s for GCs when calculating Trtot— hb)- 
For colour method, panel (a) shows the comparison between 
r (tot-HB) and r. And panels (b) and (c) represent the level 
of comparisons between a case of (tl, t) and a case of 
(tl, T( tot _HB))- We give the range for r( to t-HB)-r = 3.0 
by dashed line in panel (a). We can see all of GCs with blue 
HBR are with T( to t-HB) > r i and part of GCs with red and 



unknown HBR are with T( to t-HB) > T - On the whole, the in- 
fluence of HB stars on age determination is larger for GCs 
with blue HBR than those with red and unknown HBR. 
From panels (b) and (c), we can see that most of GCs are 
on the t 5Ctl side in panel (b), while the GCs are on both 
side of the dotted line (1:1 relation, T( tot _ H B)— Tl) in panel 
(c), and the dispersion between TYtot-HB) an d tl is smaller 
than that between r and tl. So the consistence between 
7"(tot-HB) and tl is better than that between r and tl ex- 
cepting four GCs with small TVtot— hb) an d T. These four GCs 
are with r = 5.0, and they are the GCs (NGC 6342, 6338, 
6441 and 6652) said in Section 4.2. From these two panels we 
can know that the HB stars can change the age of these four 
GCs about 1.0—2.0 Gyr and can not change age largely. Be- 
cause we fix [Fe/H]c on calculating 7"(tot-HB) and two (NGC 
6342 and 6652) of these four GCs have relative large [Fe/H]c, 
these two GCs can be affected by the age-metallicity degen- 
eracy. And the result would be T( to t-HB)3 > ' 1 " if [Fe/H]c is 
not fixed. However, for the other two GCs, they have a bit 
larger [Fe/H]c, excepting affected by age-metallicity degen- 
eracy, they also can be affected by their positions on Galactic 
plane for they being close to Galactic bulge. 

For spectrum-fitting method, just as colour method, 
panels (d), (e) and (f) show the comparisons between 
T(tot-HB), T and tl- And we give the range for r {tot _ HB )-r = 
3.0 by dashed line in panel (d). We can see all of GCs with 
blue HBR are with T( tot _ HB ) >t, and part of GCs with red 
and unknown HBR are with T( tot _ H B)>T. On the whole, the 
influence of HB stars on age determination is larger for GCs 
with blue HBR than those with red and unknown HBR. 
From panels (e) and (f), we also can see that the dispersion 
between TYtot-HB) and tl is smaller than that between r and 
tl and the GCs is more close to the dotted line (1:1 relation) 
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Figure 8. The influence of HB stars on the age determinations for colour (top panels) and spectrum-fitting (bottom panels) methods. 
Red open squares, blue open circles and green open pentacles stand for GCs with red, blue and unknown HBR, respectively. The dotted 
lines present 1:1 relation, the dashed lines indicate the range for t (tot— HB) — T = 3-0. 




Figure 9. The influence of HB stars on the age determinations for Lick-index method. Lines are the same as in Fig. 3. The solid symbols 
are the GCs with HB stars and the open symbols are the GCs without HB stars. The squares, circles and pentacles stand for GCs with 
red, blue and unknown HBR, respectively. The arrows indicate how the indices for GCs change when HB stars are removed, and we do 
not show the arrow for GCs with red HBR because the influence of HB stars for this kind of GC is very small. 
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in panel (e), so the consistence between TYtot-HB) an d Tl is 
better than that between r and tl. 

In the 10th and 12th columns of Table A2, we give the 
change of age (Arg B and Ar| B , T (tot _ HB ) — t) influenced by 
HB stars for these two methods. On the whole, the existence 
of HB stars can make age be small about 0.0—3.0 Gyr for 
this two methods. 

For the Lick-index method, the influence is displayed in 
Fig. 9, in which the open and solid symbols stand for GCs 
with and without HB stars, respectively; squares, circles and 
pentacles display GCs with red, blue and unknown HBR, re- 
spectively. The positions of GCs with and without HB stars 
change a lot for those with blue and unknown HBR and 
small for those with red HBR. For clarity, we use arrows to 
indicate the moving directions when HB stars removed for 
GCs with blue and unknown HBR whose positions change 
large, and do not use arrows to indicate the moving direc- 
tions for GCs with red HBR whose positions change small. 
From these arrows, we can see that the positions for those 
with blue and unknown HBR are shifted toward older age 
(>15.0 Gyr). All these show that the ages of GCs with blue 
and unknown HBR can be affected strongly by HB stars, 
and the ages of those with red HBR are affected slightly 
by HB stars. In the 11th column of Table A2, we list the 
change of age (Ar^ B , T( to t-HB) — t) affected by HB stars 
for Lick-index method. For the GCs with and without HB 
stars all lying inside Lick indices grids, we give the numeri- 
cal change of age. For those with HB stars lying inside grids 
and lying beyond grids when HB stars removed, we assume 
At ^ 5.0 Gyr to show the change of age even though this 
may be just the minimal change of age. For those with and 
without HB stars all lying beyond grids, we use 'uncom' in- 
stead of numerical value to express the change of age, but 
this does not mean that the age of GCs changes small. On 
the whole, the lower limit of maximal change of age is 6.0 Gyr 
due to the influence of HB stars, and the influence of HB 
stars is stronger for GCs with blue and unknown HBR than 
those with red HBR. 

5.2 Blue straggler stars 

BSs are identified as blue, luminous extensions of main- 
sequence (MS) stars (Sandage 1953), and they have been 
widely observed in stellar systems (Johnson et al. 1999; Al- 
caino et al. 2003; Piotto et al. 2004; Ahumada & Lapasset 
2007). Just as the blue HB stars, the existence of BSs can 
also mimic the presence of younger SPs (Lee et al. 2000; 
Schiavon et al. 2004; Cenarro et al. 2008). They can en- 
hance the integrated spectrum in ultraviolet and U bands 
of the GC. This phenomenon would lead to the GC being 
predicted as young or lower metallicity based on the EPS 
models without considering the BSs. Cenarro et al. (2008) 
found that higher BSs ratios can lead to smaller apparent 
ages. Xin et al. (2008) studied the BSs' influence on the in- 
tegrated properties of Large Megellanic Cloud star cluster 
ESO 121 — SC03, and found that the best-fitting values of 
age and metallicity are significantly underestimated com- 
pared to the true cluster parameters. 

Similar to the procedure of investigating the influ- 
ence of HB stars on age determinations in Section 5.2, we 
also calculate the influence of BSs on age determinations. 
The first step, we select the BSs from the CMDs of P02 



through directly visual inspection of the CMD of each clus- 
ter (Fig. 6, just as the method describing in Moretti, An- 
geli & Piotto 2008). And in this work we select 9 GCs 
(in common with our sample) with distinct BSs to study 
the quantitative influence of BSs. The parameter of T c g 
for each BSs can be estimated by the empirical re lation 
log T cff = -0.38(b BS - vbs) +3.99 (jFerraro et alj|2006l ). Fcr- 
raro et al. (2006) also indicated that the gravity of BSs was 
in the range of 4.3^ logg ^4.8. In our procedure, we choose 
log<? = 4.5 for BSs. The same metallicity of GC is adopted 
for BSs of a certain GC. Based on these parameters we can 
obtain the colours, Lick indices and spectra of BSs from stel- 
lar spectral library. Because the Bluered library is at high 
resolution, it does not provide colours. We adopt the colours 
from BaSeL library and spectra from Bluered library. 

Using the methods described in Section 3, we can study 
the influence of BSs on age determinations. Similar to Fig. 8, 
Fig. 10 shows the influence of BSs on age determinations for 
colour and spectrum-fitting methods, and T( tot _ B g s ) repre- 
sents the age of GC without BSs. For colour method, panles 
(a), (b) and (c) show the comparisons between T( tot _ B s s ), 
t and tl. We give the range for T( tot _ B g s ) — r = 2.0 by 
dashed line in panel (a). We can see that most GCs are 
with T( tot _BSs)>T, and this indicates that the existence of 
BSs can make GCs look younger. From panels (b) and (c), 
we can see the GCs lie on both side of the dotted line in 
panel (c) and lie on one side of the dotted line in panel (b) , 
so the consistence between T( tot _ B s s ) an d tl is better than 
that between r and tl. 

For spectrum-fitting method, similar to colour method, 
panels (d), (e) and (f) show the comparisons be- 
tween 77tot-BSs)! T and tl. And we give the range for 
T (tot- B Ss) —t = 2.0 by dashed line in panel (d). Just as 
the colour method, most GCs are with T( tot _ B s s ) >t. From 
panels (e) and (f), we can see that the dispersion between 
7"(tot-BSs) an d t l is smaller than that between r and tl, and 
the GCs is more close to the dotted line (1:1 relation) in 
panel (e), so the consistence between T( tot _ B g s ) and tl is 
better than that between r and tl. 

In the 13th and 15th columns of Table A2, we give the 
change of age (Arg Sa and At| Sb , T( tot _ B sa) - t) influenced 
by BSs for these two methods. The change of age is about 
0.0—2.0 Gyr for the influence of BSs, and this approximates 
to influence (0.0—3.0 Gyr) of HB stars on age determinations 
for these two methods. In theory, the T e g of HB stars is 
higher than BSs, and the existence of HB stars can make 
the spectra and colours bluer than that of BSs. But these 
two type stars have great contribution in blue and ultraviolet 
bands. In this work, we select U BVI colours and the spectra 
in the 3700-5700 A range, so the difference of influences 
by these two type stars is small in such range of colours 
and spectra, and the change of age affected by these two 
type stars is nearly the same. The difference of influences on 
age by these two type stars may be obvious for Lick-index 
method. 

Fig. 11 displays the effect of BSs on age determination 
for Lick-index method. Open triangles and pentacles stand 
for GCs with and without BSs, respectively. From these dif- 
ferent symbols, we can see that the Balmer indices of GCs 
with BSs are stronger than those without BSs, and those 
without BSs go down on the Lick indices planes. Similar 
to HB stars, the change of age (Arg Sa , i~( tot _ B s s ) — r) af- 
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Figure 10. The same as Fig. 8, but for the influence of BSs on the age determinations. 
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Figure 11. The influence of BSs on the age determinations of GCs for Lick-index method. Lines are the same as in Fig. 9. Open triangles 
are the GCs with BSs, and open pentacles are the GCs without BSs. 



fected by BSs for Lick-index method is listed in the 14th 
column of Table A2. On the whole, the lower limit of max- 
imal change of age is 5.0 Gyr due to the existence of BSs, 
and this approximates to the lower limit of maximal change 
of age affected by HB stars (6.0 Gyr). We can not say that 
the influences of this two type stars are similar, because we 



do not give the numerical value of age change for those lying 
beyond grids. From Figs. 9 and 11, we can see that HB stars 
can make H/3 stronger about 1.5A, but BSs can make H/3 
stronger about 0.5A, this demonstrate that the influence on 
age determination for HB stars is larger than that of BSs. 
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Figure 12. The influence of binary stars on age determinations 
for colour (left panel) and spectrum-fitting (right panel) methods. 
The dotted and dashed lines have the same meaning with Fig. 8. 



And the difference of influence on age by these two type 
stars is obvious for this method. 



5.3 Binary stars 

Binary stars are very common in star clusters and galaxies. 
The evolution of binary stars is very different from single 
stars, and binary interactions can also create some impor- 
tant objects and phenomenons, such as BSs (e.g. Pols & 
Marinus 1994; Tian et al. 2006), subdwarf B stars (Han et 
al. 2002, 2003), and the Gravitational wave radiation sources 
(Liu 2009) . Therefore, binary stars have the potential to play 
an important role in determining the overall appearance of 
any realistic SP. Zhang et al. (2004, 2005a) have investigated 
the effects of binary interactions on the integrated colours, 
spectra and Lick indies. And they found that the inclusion 
of binary interactions made the integrated U — B and B — V 
colours and spectra of populations bluer for various metal- 
licities, and made the H/3 greater than that without binary 
interactions. 

In this work, we investigate the effect of binary stars 
on the age determinations of GCs based on the models of 
SSPs (Zhang et al. 2004) and BSPs (Zhang et al. 2005a). 
For these two models, we adopt the Bluered library instead 
of BaSeL for spectrum-fitting and Lick-index methods, and 
BaSeL library for colour method. More description of these 
two models is given in Section 2.2. In this part we obtain 
the ages of GCs from three methods based on SSP and BSP 
models. In Fig. 12 we display the ages determined from SSP 
and BSP models for colour and spectrum-fitting methods. 
The ages of GCs obtained from BSP model are lager than 
those obtained from SSP model. Similar to the HB and BS 
stars, the existence of binary star can make the GCs look 
younger (about 0.0 — 3.0 Gyr) for these two methods. Similar 
to above, binary stars also have great contributions in blue 
and ultraviolet bands. So the influence of binary stars is 
small in the colours and spectra that we used. Fig. 13 gives 
the results for Lick-index method, in which the solid and 
dashed lines stand for BSP and SSP models, respectively. 
From them we can see that the grids for BSP model go 
up in these five panels, and the five Balmer indices of BSP 
model are greater by ~0.15 A than those of SSP model. 
Using these two models to study the ages for GCs, we can 



see that the lower limit of maximal change of age is 3.0 Gyr 
for the existence of binary stars. 



5.4 a enhancement 

The existence of a— enhancement can make stars on the 
isochrones hotter throughout all evolutionary phases, and 
the effect increases w ith metallicity. Galac t ic GCs are found 
to b e a— enhanced (|Thomas et al.l 120031 ; iMaraston et al. 1 
120031 ; iLee fc Worthevl 120051 ). Just as the other factors, 
the existence of a— enhancement also can mimic young 
SPs and make Galactic GCs look younger. For the Lick- 
index method, as said in Section 3.2, the Balmer indices 
would be affected by a— enhancement, so the age deter- 
mined by Lick-index me t hod c an be affected by th is factor . 
IMaraston fc Strombackl (|201ll ) and iThomas etaT] (|2011bl ) 
have selected some Lick indices and used the x 2 ~ fitting 
technique to derive the parameters for GCs, and they 
found that the derived ages agreed well with those ob- 
tained by CMDs when takin g account the a— enhancement. 
I Thomas et al.l (|2003l . l2011al ) have presented the whole set 
of Lick indices of SSP models with variable a/Fe ratio, this 
allows us to study the influence of a— enhancement on age 

determination for Lick-index method. 

We use the Lick indices of lThomas et al.l l|2011al ) models 
to study the influence of a— enhancement on age determi- 
nation for Lick- index method, and two types of a/Fe ratio 
([a/Fe]= 0.0 and [a/Fe] = +0.5) are chosen. Fig. 14 gives 
the result, from them we can see that a— enhancement can 
make the Balmer indices stronger, especial for high metal- 
licity, and has little effect on [MgFe]' which consists with 
above said (Section 3.2). Therefore the a— enhancement can 
affect the age determination and has little effect on metal- 
licity determination for Lick-index method. We also can 
see that the influence of a— enhancement on H/3 and H7F 
is small at the range of [Fe/H] < —0.7. In this work, the 
ages for GCs are computed as the weighted mean of the 
ages derived from five Balmer indices versus [MgFe]' grids 
and most of the GCs are in the range of [Fe/H] < —0.7. 
Therefore, the a— enhancement has little effect on age de- 
termination of GCs for Lick-index method. But for the 
study of high-metallicity SP systems, the a— enhancement 
should be taken into account, because the influences of 
q— enhancement on these five Balmer indices are all signifi- 
cant for high-metallicity range. 



6 SUMMARY 

We have investigated the utility of three methods (colour, 
Lick-index and full spectrum-fitting) with the same EPS 
models (BC03 models) for estimating ages and metallicities 
of Galactic GCs. We also compared our results with those 
estimated in the literatures from other methods, and the 
main results of this study are as follows. 

(1) Our results show that Galactic GCs are almost old 
metal-poor SP systems. 

(2) Metallicities determined from these three meth- 
ods are in agreement with those of literatures in the en- 
tire range spanned by GCs. The Lick-index method is suit- 
able to obtain metallicity for SP systems in the range of 
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Figure 13. The influence of binary interactions on the age determinations on Balmer vs. [MgFe] panels for some GCs. Solid and dashed 
lines stand for BSP and SSP models, respectively. The dots stand for GCs. 




Figure 14. The distribution of GCs on the Balmer indices and [MgFe] planes based onflhomas ct al. (2011a]) model. The lines are the 
same as in Fig. 2, the gray solid and red dashed lines stand for [o/Fe]= 0.0 and 0.5, respectively. 
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— 1.5 <[Fe/H]< —0.7 and spectrum-fitting method is suit- 
able to study metallicity for SP systems in the range of 
-2.3 <[Fe/H]< -1.5. 

(3) There exist some discrepancies between our ages 
and the literatures, we can not directly say which method 
is more suitable for age determination. Through comparing 
with the literatures we find our results are younger than the 
literatures on average, especially for GCs with blue HBR. 
Many factors can affect the age determinations and make 
GCs become younger (such as the HB morphology, BSs, bi- 
nary stars, a— enhancement) . 

(4) We use Vazdekis and Maraston models to investi- 
gate whether our results are dependent on EPS models. The 
comparisons of derived metallicities and ages for individual 
GCs show that three models make different parameter pre- 
dictions. However, the whole tendency for these three meth- 
ods of these three EPS models is nearly the same. All these 
indicate that our above conclusions are independent of the 
EPS models. 

(5) In this work we use the old Galactic GCs to test 
these three methods, and our results may hold for old SP 
system. 

(6) We also study the quantitative influences of HB, 
BSs and binary stars on age determinations for three meth- 
ods. The existence of all these stars can make the GCs look 
younger. For the colour and spectrum-fitting methods, the 
age can be underestimated about 0.0—3.0 Gyr, 0.0—2.0 Gyr 
and 0.0— 3.0 Gyr due to influences of HB, BSs and binary 
stars, respectively. And for Lick-index method, the lower 
limit of maximal change of age is 6.0 Gyr, 5.0 Gyr and 
3.0 Gyr due to influences of horizontal branch, blue strag- 
gler and binary stars, respectively. 

(7) For the Lick-index method, we also investigate the 
influence of a— enhancement on age determination, and find 
that the a— enhancement has little effect on age determina- 
tion for Galactic GCs. 
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APPENDIX A: THE PARAMETERS OF 
GALACTIC GCS AND SOME INFLUENCE ON 
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Table Al. The main basic parameters for GCs. The first column gives the names for Galactic GCs. The 2nd and 3rd columns are 
right ascension (RA) and declination (DEC) (epoch J2000). Galactic longitude (L) and latitude (B) in degrees are listed in the 4th and 
5th columns. Distance from Sun and Galactic center are shown in the 6th and 7th columns. The 8th column gives the reddening value, 
E(B — V). The parameter describing the HB morphology is listed in the 9th column. All of these data are adopted from Harris' catalogue. 
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